NASA MEMO 3-9-59A 


— NASA MEMO 3-9-59^ 

:/ : 

NASA 

MEMORANDUM 


A WIND-TUNNEL INVESTIGATION OF THE STALL-FLUTTER 
CHARACTERISTICS OF A SUPERSONIC-TYPE 
PROPELLER AT POSITIVE AND 
NEGATIVE THRUST 

By Vernon L. Rogallo and Paul F. Yaggy 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 

WASHINGTON 

May 1959 




NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


MEMORANDUM 3-9-59A 


A WIND-TUNNEL INVESTIGATION OP THE STALL-FLUTTER 
CHARACTERISTICS OF A SUPERSONIC-TYPE 
PROPELLER AT POSITIVE AND 
NEGATIVE THRUST 

By Vernon L. Rogallo and Paul F. Yaggy 


SUMMARY 


An investigation has been made utilizing a three-blade , 10 -foot- 
diameter, supersonic- type propeller to determine propeller flutter charac- 
teristics. The particular flutter characteristics of interest were (l) the 
effect of stall flutter on a propeller operating in positive and negative 
thrust, (2) the effect of stall flutter on a propeller operating with the 
thrust axis inclined, and (3) the variation of vibratory blade shear 
stresses as the stall flutter boundary is penetrated and exceeded. Thrust 
and power measurements were made for all test conditions. Wake and inflow 
surveys were made, when appropriate, to define the thrust and torque dis- 
tributions and the magnitude of the inflow velocity. Stress measurements 
were made simultaneously to obtain the propeller flutter and bending 
response. 


It was found when operating both in the positive and negative thrust 
regions that, for most cases after the onset of flutter, the magnitude of 
flutter stresses at first increased rapidly with section blade angle, 

3, after which further increases in p resulted in only a moderate 
increase or a reduction in stress. Thrust-axis inclination up to the 
limit of the tests (angle of attack of 15° and dynamic pressure of 4o psf ) 
appeared to have no effect on stall flutter. The stall flutter stresses 
were found to be directly associated with the section thrust characteris- 
tics of the blades. The onset of flutter was found to occur simultaneously 
Ihe divergence of the section thrust variation with blade angle from 
linearity for stations outboard of the blade 0 . 8 -radius station. The 
maximum flutter stresses appeared to be a function of the maximum section 
thrust obtained at or in the vicinity of the blade 0 o 8- radius station. 

In an attempt to correlate two-dimensional airfoil data with three- 
dimensional data to predict the stall angle of attack (divergence of the 
section thrust) of the blade sections, it was found that no consistent 
correlation could be obtained. Also, a knowledge of the inflow conditions 
appeared to be insufficient to account for differences in airfoil charac- 
teristics between the two-dimensional and the three-dimensional cases. 
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INTRODUCTION 


Propellers suitable for high subsonic and transonic speed aircraft 
have blades with very thin airfoil sections. Blades which have thin 
airfoil sections are flexible both in bending and in torsion. To absorb 
the high power available from turboprop engines during take-off and early 
climb , these thin, flexible blades must often operate at blade angles 
above those at which blade-section stall v.ould occur. The dynamic char- 
acteristics of the blades coupled with the required operating conditions 
is conducive to a blade flutter phenomenon in torsion termed stall flutter. 
Stall flutter has been so violent, in some cases, that the design was 
precluded by structural limitations. 

Semiempirical methods for prediction of the stall flutter boundary 
(i.e., the operating conditions at which stall flutter will occur) have 
been developed. However, the general applicability of empirical methods 
is not evident unless they are based on ar adequate understanding of the 
characteristics and mechanics of stall flutter. Some characteristics yet 
unknown and believed worthy of investigation are as follows: 

1. The variation of vibratory blade stresses as the stall flutter 

boundary is penetrated and exceeded. 

2. The effects of stall flutter on propeller operating charac- 

teristics . 

3. The effects of blade twist, power input, design lift coefficient, 

solidity, number of blades, etc. , on the stall flutter boundary. 

1+. The applicability of the resell ts of scale-model tests to full- 

scale design. 

To provide information for items 1 and 2, full-scale tests of a 
supersonic-type propeller have been made. The tests were made in the 
Ames ho- by 80-foot wind tunnel which, because of its large size relative 
to the size of the propeller, minimized tunnel-wall interference. The 
propeller was selected because its dynamic characteristics were such that 
the blades were susceptible to stall flutter over a wide range of blade 
angles and rotational speeds. Tests were made with the propeller operat- 
ing at both positive and negative thrust with the thrust axis parallel to 
the wind stream. Tunnel airspeeds corresponded to a range from near zero 
through early climb or landing and rollout speeds. Tests were made also 
at positive thrust conditions with the thrust axis inclined to the wind 
stream for early climb speeds. Vibratory shear stress characteristics of 
the blades and thrust and power character:. sties of the propeller before 
and during stall flutter were measured for all operating conditions. 

Inflow velocity and wake surveys were made for positive thrust conditions. 
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A correlation has teen made of the onset of flutter as well as the maximum 
flutter stresses with measured aerodynamic data. In addition, an analysis 
has been made to determine the applicability of two-dimensional airfoil 
section data to predict the three-dimensional stall angle , defined as the 
divergence of the section thrust variation with blade angle from linearity 


NOTATION 


b blade- section chords in. 
c^ section drag coefficient 
c^ section lift coefficient 

c^ two-dimensional lift-curve slope, per deg 
Cj. elemental blade-section thrust coefficient. 




elemental blade-section torque coefficient, 

P 

power coefficient, — «— 

Pn 3 D 5 



torque coefficient. 


Q 

Pn 2 D 5 


T 

C T propulsive thrust coefficient, — - — 

pn 2 D 4 

D propeller diameter, ft 

V 

J propeller advance ratio, 

Mg blade- section Mach number 

N propeller rotational speed, rpm 

n propeller rotational speed, rps 

P power input, ft-lb/sec 

Q propeller torque, ft- lb 

r radius to any blade section, in. 



k 


R tip radius, in. 

T propulsive thrust, lb 

t maximum blade-section thickness, ir . 

V velocity of free-stream tunnel air stream, fps 

oo 

V inflow velocity at propeller plane induced by blade loading, fps 
total inflow velocity at propeller plane, Vj + cos Oq, fps 

r 

X fraction of blade radius, ^ 

blade-section angle of attack, deg 

maximum variation (first-order) of a, 

geometric angle of attack of the thrust axis, deg 

section blade angle, measured from a plane perpendicular to the 
thrust axis, deg 

section blade angle at the 0 . 7“ radius station 

Vi 

geometric angle of advance, tan ’ tn p x > ^- e 6 
mass density of air, slugs /cu ft 

MODEL AND APPARATUS 
Dynamometer and Propeller 



The propeller and test rig used for the investigation are shown in 
figure 1 as they appeared mounted in the 4 > by 80 -foot wind tunnel. The 
propeller was driven by two 1500 horsepower (3000 rpm) frequency- 
controlled electric motors coupled in tandem and connected directly to 
the propeller. The three-blade, 10 -foot-d Lame ter propeller used for the 
investigation was designated as Curtiss C636D-A3X. The blades, designated 
design number 109640 , were solid dural and had symmetrical l6-series air- 
foil sections. The geometric characteristics of the blades are given in 
figure 2 . This type of propeller has been termed "a supersonic propeller" 
since at the design cruise speed the blade airfoil sections are operating 
at supersonic speeds. 

The propeller blade angle was controlled by remote operation of the 
standard integral pitch control mechanism contained in the C636 hub. 
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Instrumentation 


The blades were instrumented with Baldwin- Southwark type SR- 4- CB10 
strain gages as shown in figure 3(a). The gage locations are shown in 
figure 3(b). All blades had bending gages at 20 inches from the thrust 
axis and shear (torsion) gages at b-2 inches from the axis. These positions 
were those where the maximum vibratory stress had been calculated to occur 
and the readings from these gages provided the phase relation and relative 
stress level of the blades to each other. The master blade had additional 
bending gages at 30 , 40, and 50 inches from the thrust axis to define the 
radial variation of the flatwise bending and an additional shear gage at 
24 inches from the axis to identify the torsional mode. The output signals 
of the gages, appropriately attenuated, were simultaneously directly 
recorded on a Consolidated 36 -channel, moving-coil- type oscillograph using 
type 7-218 galvanometer elements. Only the vibratory stresses were 
recorded; the steady stresses were eliminated by blocking condensers. 

Propeller blade angle was measured by means of a potentiometer 
attached to the blade index head of the integral pitch-control mechanism 
of the propeller. The potentiometer output was recorded on a direct- 

reading Brown potentiometer modified and calibrated to indicate the blade 
angle in degrees. 


The propulsive thrust of the model was measured by the longitudinal 
force scales of the wind-tunnel balance system* 

Measurements of the flow field were made in both the wake of the 
propeller and just upstream of the propeller. The wake surveys were made 
y means of an eight- tube directional pitot-static rake located 19 - 3/8 
inches aft of the propeller plane. The rake is shown in figure 1. From 
e readings of this rake, both propeller thrust and torque distributions 
were obtained. The inflow surveys were made by means of a rake composed 
of three static and three total- head probes. The probes were located 
0 - 3 /I 0 inches forward of the propeller plane and at the 0 38 - 0 58 - 

and 0.78- radius stations. This rake is also shown in figure 1. *Both 

rakes were connected to multiple- tube manometers which were recorded 
phot ographi cally . 


The power input to the electric motors was measured by standard 
polyphase, laboratory- type wattmeters. The wattmeters were calibrated 
in terms of propeller torque by suspending the test rig in a cradle which 
allowed torque measurement. 


The rotational 
Berkeley electronic 
in the motors . The 
by means of a pulse 
magnet, attached to 
rig. 


speed of the propeller was observed by means of a 
counter connected to a tachometer generator located 
speed was also recorded on the oscillograph record 
which was generated once per revolution by an alnico 
the propeller hub, passing a fixed coil on the test 
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TESTS 

Propeller Nonrotating Vibrational Characteristics 


The fundamental torsional frequencies of the nonrotating blades were 
determined with the propeller installed on the test rig. The blades were 
excited by electromagnets. The force and frequency of oscillation of the 
magnetic fields were controlled by an elecironic power supply incorporating 
an audio- oscillator. The frequency of eacf blade was measured with the 
excitation applied at that blade. The tests were made at a blade angle 
of 13 5° measured at the 0.7" ra< 3-ius station. The fundamental torsional 
frequencies were found to be 113.5, U2.5, and 113-2 cycles per second 
for blades 1, 2, and 3, respectively. 

The decay technique was used to determine the damping characteristic 
of each blade (i.e., the excitation power was cut off abruptly and the 
decay of the amplitude of vibration was recorded on the oscillograph) . 

From these records for a given amplitude o:‘ excitation, the logarithmic 
decrements, 5, were determined and found to be 0.0055, 0.0050, and 0.0051 
for blades 1, 2, and 3, respectively. 


Wind-Tunnel Tes :s 


The range of test conditions is tabulated in table I. Because of 
model power and blade stress limitations, the propeller could not be 
operated over the entire propeller rotational speed and blade angle ranges 
for all tunnel airspeeds. The test conditions generally began just prior 
to the onset of stall flutter and the stall flutter boundary was penetrated 
to the maximum conditions permitted by stress and power limitations. A 
few runs were made also to determine the propeller operating characteris- 
tics beginning at a condition far removed from flutter and proceeding 
through the flutter boundary to the maximum blade stress permissible. 

Two techniques were used to penetrate the stall flutter boundary . 

The first was to establish a constant propeller rotational speed and to 
vary the blade angle while a constant tunr el velocity was maintained. 

The second was to establish a constant blEde angle and to vary propeller 
rotational speed while a constant tunnel velocity was maintained. The 
first technique was found to be the most e atisfactory from an operational 
standpoint and was employed throughout moE t of the tests. 

The data were all recorded simultaneously by photographing instrument 
panels and multitube manometers with cameras which were synchronized with 
the recording oscillograph. 
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REDUCTION OF DATA 


The effect of blade instrumentation on the airloads and resultant 
stresses is believed to be small, since the results of previous investi- 
gations on full-scale propellers have shown that the agreement in magnitude 
of blade stresses in a fully instrumented blade and blades of limited 
instrumentation were within the accuracy of the measurements. Tests were 
made with and without the inflow survey rake, under as nearly identical 
conditions as practical, to determine the effect of the rake on the onset 
of flutter and on the magnitude of the flutter stresses. The rake appeared 
to have no effect on either the onset of flutter or the magnitude of the 
flutter stresses. No corrections for tunnel-wall interference have been 
made since the effects of tunnel-wall interference are believed to be 
small because of the large ratio of the tunnel cross-sectional area to 
the disk area of the test propeller. 


Oscillograph Records 


Since the complete strain-gage system had been calibrated in terms 
of stress per unit deflection of the recording galvanometers, the blade 
vibratory stresses could be obtained directly from the oscillograph 
records. One record was taken for each test condition. The average length 
of record was 60 inches and was recorded at a speed of 4o inches per 
second. The value of stress at any gage location for a given test con- 
dition was taken to be the maximum stress indicated on the oscillograph 
record for that condition. The frequency of the stress variation was 
obtained directly from the oscillograph records by means of timing lines 
spaced at intervals corresponding to 0.01 second. 


Propeller Wake-Survey Data 


Blade element thrust and torque coefficients were computed by the 
method given in reference 1. This method required measurements of total- 
head pressure and flow angle immediately aft of the propeller plane in 
the absence of the propeller and in the presence of the propeller. 

Previous investigations (ref. 1 and others) have shown this to be an 
acceptable means of measuring thrust and torque (it is realized that some 
error exists when the blade sections are operating at a stalled condition). 


Propeller Inflow Data 


The inflow velocity at the propeller disk was determined by linear 
interpolation between the value of velocity measured at the inflow rake 
and value of velocity measured at the rake in the propeller wake. This 
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procedure is believed to be in accordance with momentum theory. Both 
survey rakes were calibrated prior to the test. The data obtained at 
low forward airspeeds were found to be inaccurate, in particular at the 
outboard blade stations, because of high flow angles caused by the con- 
traction of the inflow stream (the inflow survey tubes were not of the 
directional type ) . Thus , the data presents d herein are for a forward 
speed range of 58.5 to 183 *5 fps. 


Force Test Dat^ 


The propulsive thrust is defined as the longitudinal force with the 
propeller operating less the longitudinal force with the propeller removed, 
both at the same tunnel air-stream velocity. 


RESULTS AND DISCUSSION 


The results of this investigation and discussion pertinent thereto 
will be divided into three main sections, (l) typical flutter stress 
characteristics at both positive and negative blade angles (included 
herein are presentations and/or discussions of blade phasing at flutter, 
flutter boundaries, wake-excited flutter, and the effects of thrust-axis 
inclination on flutter), (2) section thrust and inflow velocity character- 
istics as determined from survey rake measurements (included herein are 
correlations of the aerodynamic data with flutter stress data, and an 
analysis to determine the applicability of two-dimensional airfoil section 
data to the prediction of the three-dimensional stall angle), and (3) pro- 
peller performance characteristics as determined from the tunnel force 
scales and power measurements (these include absolute values of both 
thrust and power for a typical case at positive and negative thrust, in 
addition to the thrust and power coefficients for all test conditions of 
the investigation) . 


Flutter Stress Characteristics 


Flutter at positive thrust.- The flutter stresses as presented herein 
apply to all of the blades for' a given test condition, since, for all con- 
ditions of sustained flutter, all three blades exhibited very nearly 
identical flutter characteristics. This could be expected from the results 
of the nonrotating tests which indicated '-hat the blades were dynamically 
similar. Typical vibratory shear stress characteristics for the case of 
positive blade angle operation are shown .n figure 4 for near- zero tunnel 
velocity and at moderate tunnel velocities. Aside from the obvious effects 
of increasing tunnel velocity (i.e., causing the initial shear stress rise 
to occur at higher blade angles), the stress rise characteristics for all 
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tunnel velocities are similar for a given rotational speed. Note that 
for most constant rotational speed runs with the exception of those for 
which stress peaks could not he reached (limitations in power, allowable 
stress, etc.) the stress at first increased rapidly with 3, after which 
further increases in 3 resulted in only a moderate increase in stress 
and in some cases a reduction in stress. It was noted during the test 
that at the first indication of vibratory shear stress there was no con- 
sistency in the response of the three blades. It often occurred that one 
blade was experiencing a low value of vibratory shear stress (less than 
±1000 psi) when a second would suddenly trigger and experience a somewhat 
higher value of vibratory shear stress. The first indication of vibratory 
shear stress for all blades was depicted by a wave form sinusoidal in shape 
but having an apparent beat frequency. This pattern was neither of a 
constant amplitude nor of a constant beat frequency. As the control 
parameter (rpm or 3) was further increased, the beating tended to dis- 
appear when the value of shear stress increased to the order of ±1000 psi. 
In general, this occurred at about the same time for all blades. Further 
increase in the control parameter resulted in higher levels of shear stress 
of nearly equal magnitude for all blades and the wave form of the oscilla- 
tion was of a sustained amplitude. 

It is of interest to note the marked difference in the stress rise 
characteristics for propeller speeds below 1200 rpm, in contrast to those 
at and above this value. A consideration of two-dimensional flow indicates 
that these stress characteristics might be expected. Although the 
mechanics of flutter are not precisely known, it is believed that the 
variation of the lift forces on an oscillating airfoil is dependent upon 
the relative length of the repeating vortex pattern in the wake of the 
airfoil as compared to the chord length (i.e., the number of chords that 
the vortex pattern moves downstream in one cycle of vibration defines, in 
part, the forces on the airfoil). A parameter describing the geometry of 
the flow field is commonly referred to as the "reduced frequency parameter" 
and is defined by 


where 

id frequency of the vibration, radians/sec 

b semi chord, ft 

V local air velocity, ft/sec 


Investigations have shown that the critical mean angle of attack at which 
energy is absorbed from the air stream is a function of this parameter. 

It has been found convenient to use the reciprocal of the reduced fre- 
quency parameter to denote a flutter speed coefficient. A, written as 
(ref. 2) 


A = 


V 


bid, 


a 
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Experimental results indicate that fluttei will not occur, regardless of 
the airfoil angle of attack, if the value of A is less than unity. For 
the case of the propeller this coefficient varies in value along the blade 
because of the differences in local air velocity. The value of A for 
the rotational speeds below 1200 rpm (fig. 4) was generally less than 
unity over the major portion of the blade. 

Flutter frequency and blade phasing .- The torsional frequencies of 
the blades during flutter were found to be near the values of the torsional 
natural frequencies obtained during the ncnrotating tests. However, since 
all blades differed in frequency with resject to each other (maximum 
variation of frequency was from 109 to ll4 cps) and drifting occurred for 
each blade, no phase relations could be established. It was noted that a 
phase shift as great as l 80 ° often occurred within a period of three- 
tenths of a second. These results are not in agreement with those reported 
in reference 3 > wherein definite phase relations were observed to occur. 

Flutter boundaries .- The flutter boundaries based on a constant shear 
stress level of ilOOO psi , established from the data of figure 4 and other 
similar data of this investigation, are shown in figure 5- Included is the 
static flutter boundary established from the data obtained at the WADC 
Laboratory for the test propeller. The choice of a criterion by which to 
establish flutter boundaries is somewhat arbitrary. The constant stress 
level criterion was selected, based to a great extent upon the character- 
istics of the onset of flutter which were discussed previously. For the 
test propeller, a value of ±1000 psi shear stress was generally indicative 
of a transition from a varying amplitude cf the shear stress to the first 
sustained value of stress encountered. Other investigators have also 
used ±1000 psi as a criterion for flutter boundaries for aluminum blades. 

Wake-excited flutter .- During the tests, flutter at low blade angles 
was observed on several occasions. Flutter occurred for only two rota- 
tional speeds, however, namely 1400 and 1830 rpm. Typical stress charac- 
teristics are shown in figures 6 and 7(a). This phenomenon has been 
experienced previously on other propellers and also on the test propeller 
at the WADC Laboratory. Flutter occurring at low blade angles has been 
termed "wake-excited flutter." Wake-excitsd flutter was characterized by 
a sustained value of stress nearly equal in magnitude for each blade and 
having a definite blade phase relationship. All blades were in phase at 
luOO rpm, and at 1400 rpm the blades were L20^ out of phase with respect 
to each other. This phasing is in agreement with the phase criterion 
for torsional flutter as established in reference 3 . 

Flutter at negative thrust .— Typical /ibratory shear stress charac- 
teristics for the reverse thrust regime ar; presented in figure 7 . For a 
given forward velocity, the stress rise as blade angle is changed begins 
at nearly the same blade angle for all propeller rotational speeds whereas 
at positive thrust the blade angle corresponding to the onset of flutter 
is more definitely a function of rotational speed. The rate of stress 
rise during negative thrust operation begins as it did in positive thrust 
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operation and in general the stress variations are similar. After the 
stresses reached their maximum value , further decreases in blade angle 
up to the maximum test limit of -17.5° were accompanied by an appreciable 
reduction in the absolute stress level for a given rotational speed. This 
trend is of significance since the implication is that operation in the 
negative thrust regime is possible with a propeller having marginal 
flutter characteristics. Thrust characteristics in this regime will be 
shown subsequently. 

Flutter at positive thrust with thrust axis inclined .- Because of 
the oscillatory variation of the blade angle of attack when the thrust 
axis is inclined, it would be expected that stall flutter would be induced 
at a lower propeller blade angle than it would for the uninclined case at 
the same operating condition. Further, results obtained on two-dimensional 
airfoil sections show that the hysteresis in the lift characteristics 
depends on magnitude and frequency of oscillations in the section angle 
of attack (ref. 4 ). This suggested that the oscillatory angle- of- attack 
variation induced by thrust-axis inclination could have a pronounced 
effect on propeller stall flutter. However, figure 8 shows representative 
results of tests with thrust-axis inclinations up to 15° and corresponding 
tunnel velocities as high as 183-5 fps which indicate little change in the 
stall flutter characteristics relative to the uninclined case. In order 
to gain an insight into why thrust- axis inclination was not influential 
in altering the flutter characteristics from the uninclined case, computa- 
tions were made to determine the Ax-g variation (extreme cx,g change in 
blade angle of attack around the disk) for a typical case. These computa- 
tions were based on strip analysis procedure for the following conditions: 
rotational speed of 1600 rpm, blade angle of 30 ° , thrust-axis inclination 
of 15 °, and a tunnel velocity of 183.5 fps. This condition was selected 
to provide an indication of the Axg variation for the test condition 
where the propeller blades were near the onset of flutter (±350 psi shear 
stress at zero thrust-axis inclination) , The results obtained indicate 
a variation of Axg from 5*65° at the propeller 0 . 3- radius station to 
0 . 75 ° at the 0 . 8 - radius station. These high values of Ax-g are confirmed 
by the relatively large magnitude of the measured once-per- revolution 
vibratory bending stresses (±5000 psi). The values of Aag are of 
significant magnitude and suggest that the blade would be well into the 
flutter region during a portion of a revolution, while diametrically 
opposite this position the blade would be appreciably removed from the 
flutter region if stall flutter phenomenon occurred over a large portion 
of a blade. Thus, an explanation as to why thrust-axis inclination does 
not influence stall flutter cannot be made without associating the aero- 
dynamic characteristics of the blade sections with the flutter character- 
istics. A study of the wake survey data to be discussed subsequently, 
for the uninclined case, revealed that for all propeller operating 
conditions investigated, only the outer portion of the blade (approxi- 
mately 20 percent) exhibited stall phenomenon. As indicated above, Ax-q 
over this region would be less than 0.75° and- of insuff icient magnitude 
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to significantly alter the stall flutter characteristics . The tendency 
of the inner portions of the blade to remain unstalled at angles of attack 
far in excess of what would be predicted from two-dimensional airfoil data 
will be discussed in the latter part of this report. 


Flow Surveys at Proximity of Propeller 


Correlation of onset of flutter with measured aerodynamic data .- In 
order to obtain a correlation of the onset of flutter with the aerodynamic 
phenomena occurring during flutter, comparisons have been made of the 
section thrust coefficient obtained from wake surveys and the blade 
torsional stresses. These comparisons have been limited to a range of 
rotational speed from 1200 to 1800 rpm which is believed to be the range 
of greatest interest for the propeller diarieter tested. This range is 
also above that wherein the effects of reduced frequency are of signifi- 
cance and below that wherein the effects of Mach number become significant. 
Typical results of the variation of the blade-section thrust coefficient 
with blade angle are shown in figure 9* If may be noted that the blade 
angle for onset of flutter, as indicated by the initial rise in shear 
stress in figure 4, very nearly corresponds to the blade angle at diver- 
gence from linearity of the thrust variation at the outboard sections for 
a given value of J (see fig. 9)- The divergence at the 0.78 an d 0.90 
stations occurs simultaneously with the onset of flutter (little or no 
indication of stall occurred for stations : nboard of the 0.64-radius 
station for all test conditions). This divergence is believed to be 
associated with flow- separation phenomena at or preceding blade-section 
stall. It appears then from the foregoing that stall flutter is dependent 
upon the aerodynamic characteristics of airfoil sections at or near the 
stall angle and that only a single outboarc. section need be considered 
in the prediction of stall flutter. 

Correlation of the maximum flutter stress with measured aerodynamic 
data .- Comparisons of the section thrust coefficient with flutter stress 
variations indicate that the maximum flutter stress is associated with 
the maximum thrust of the blade sections. The stress level appears to 
reach its maximum value when the maximum thrust is obtained at or near 
the vicinity of the 0.8-blade- radius station. This correlation can be 
made by comparison of figures 4(e) and 9 (a] and (b) for operating con- 
ditions at a moderate forward velocity and figures 4(a) and 10(a) for a 
case near zero forward velocity. For rotal ional speeds well within the 
flutter boundary extremities (i.e., 1200 to 2200 rpm) after the maximum 
thrust had been realized over the outer 20 percent of the blade, further 
increases in p generally resulted in a reduction in thrust and a nearly 
sustained or a reduced value of absolute si ear stress. 

Theoretical prediction of stall flutter .- The foregoing discussion 
intimates that the onset of stall flutter cccurs with the stalling (devia- 
tion of the section thrust variation from linearity) of the outboard 
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, _ „ whinh were computed for the same conditions. 

Ts'ZZlZTllrfA use of the blade-element theory and the^ 
two-dimensional lift curves in figure 13 (a) (primarily from ref. 5). The 
average lift-curve slope from = 0 to 6 was used. 

For the 0 38 and 0.58 stations, where wake rake data showed the 

curve slopes an? those calculated^by^use^of^nieasured^inflow^data.^On^ 

thxs basis alone, Qf at the 0.78 station could 
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Propeller Aerodynamic Performance 


No effect of flutter on propeller tlrust and power characteristics 

f s R( , can t) e seen by comparison cf figures 4 and 15- I -f 

was found, as can be see y ter stresse£ rise abruptly and obtain rela- 

noted that, although the fl . discontinuity in either 

tively high values with _ increasing p there is no discomm y 

the thrust or power vanatio . ^8-radius station were obtained 

by interpolating between values indicated by the wake rake tubes at sta- 

tions 0.50 and 0.64. . — — " 
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Because of this apparent lack of effect of flutter on thrust and 
power characteristics, it was considered plausible to present all of the 
thrust or power performance data on a single plot. The thrust and power 
coefficients are presented in figure 1 6. This figure serves to show the 
ranges of test variables covered during the tests. In addition, the 
repeatability of the data is also indicated since these data represent 
values obtained during different phases of the tests with duplications 
of given test conditions. 


CONCLUDING REMARKS 


The results of the propeller stall flutter investigation indicate 
that stall flutter can be associated with the aerodynamic characteristics 
of the blade sections. The onset of flutter was found to occur simulta- 
neously with the divergence of the section thrust variation with blade 
angle from linearity for sections outboard of the blade 0.78-radius 
station. The maximum flutter stress appeared to be a function of the 

maximum thrust obtained at or in the vicinity of the blade 0.8-radius 
station. 

In regard to the prediction of stall by use of two- dimensional air- 
foil data, it was found that there was a definite lack of correlation 
between two- and three-dimensional data. The apparent differences in 
c ^a f ? r the ^ wo cases were found to be quite significant and the differ- 
ences in section-stall angle of attack for the two cases were phenomenal. 
For example, where two-dimensional data would indicate stall to occur in 
the vicinity of 6° to 8°, the airfoil sections at the inboard sections 
of the blade were operating at angles of attack as high as 20° without 
appreciable stall. A knowledge of the inflow conditions appeared to be 
insufficient to account for differences in airfoil characteristics between 
the two- and three-dimensional cases. 

It appeared reasonable, for the test propeller, to use the constant 
stress criterion to establish flutter boundaries since the first sustained 
value of stress appeared to occur at the same stress value (±1000 psi) for 
all test conditions. The correlation of the flutter stress data with 
aerodynamic data obtained from wake surveys suggests that a flutter bound- 
ary representative of the onset of flutter could be more appropriate if 
based on thrust divergence rather than on a constant stress level or rate 
of stress rise, if such divergence could be predicted adequately. (it is 
assumed that the dynamics of the blade are such as to support flutter.) 

Since no definite phase relationship between blades at flutter 
occurred for the subject investigation, with the exception of the case 
of wake-excited flutter, the results with respect to blade phasing at 
flutter are contradictory to previous tests (ref. 3). 

Flutter appeared to have no influence on either thrust or power 
characteristics. 
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The trend found in the negative thrust region for the magnitude of 
the stress to reduce with increasing negative (3 after a maximum flutter 
stress was obtained may be of importance for propellers used as an air- 
plane braking force. Further tests should be made over an extended nega- 
tive P range to study this trend and also to determine the possibility 
of stress reduction by a rapid reduction oi P« 

Thrust-axis inclination up to 15° appeared to have no effect on stall 
■flutter. An analytical study made to determine the £cl b due to thrust- 
axis inclination together with a knowledge of the region of the blade 
stall from wake survey measurements indicated that the value of ^cc-g over 
the stalled portion of the blade was of such a low magnitude as to be 
insignificant in altering the flutter. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 8, 1958 
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TABLE I.- TEST CONDITIONS 


1 Voo, fps 

OQ, leg 

^o.yR> 

N, rpm 

Near zero 

0 

18 to 42 

550 to 2200 

58.5 

0 

18 to 42 

800 to 2000 

88.5 

0 

18 to 42 

800 to 2200 

115.5 

0 

18 to 42 

800 to 2200 

143.8 

0 

19 to 42 

800 to 2000 

183.5 

0 

-17.5 to 42 

1200 to 1800 

Near zero 

0 

-17.5 to -10 

550 to 2200 

58.5 

0 

-17.5 to -2 

600 to 1800 

88.5 

0 

- 17.5 to -8 

600 to 1800 

115.5 

0 

-17.5 to 2 

600 to 1800 

143.8 

0 

-17.5 to 2 

600 to 1800 

88. 5 a 

10 

18 to 42 

600 to 1800 

115. 5 a 

10 

22 to 42 

800 to 1800 

143. 8 a 

15 

21 to 42 

1000 to 1800 

183.5 

15 

23 to 40 

1200 to 1800 


a Data not presented for these conditions. 
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Figure 1.- The propeller and test rig installed in the Ames 40- by 80-foot wind tunnel. 



(a) Arrangement of strain gages on blade. 
Figure 3.- Instrumentation of test propeller. 
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Figure 4- # - Flutter characteristics at positive thrust; 





0>) V ro = 88.5 fps 

Figure 4.- Continued 






(c) v m = 115.5 fps 

Figure 4 .- Continued 
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Figure 4.- Continued. 
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Figure 5 .- Flutter boundaries for constant forward velocities; shear stress level criterion 

(±1000 psi); a f = 0°. 








( b ) v oo - 1^3.8 fps 
Figure 7.- Continued 
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Figure 8.- Comparison of flutter characteristics at positive thrust with 
thrust axis inclined and noninclined; = 183.5 fps. 
























(Id) Survey tube locu.tion; x 
Figure 10.- Concluded. 



















P.7R* de 6 



(k) V ro = 115.5 fps; x = 0 

Figure 11.- Continued. 
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(a) Lift characteristics. 

Figure 13.- Aerodynamic characteristics of symmetrical l6-series airfoil 

sections for several Mach numbers. 
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Figure lR.- Comparison of the computed and uake measured variation of the 
section thrust coefficient as a function of the propeller blade angle 
for several values of advance ratio; a G == 0°. 
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Figure 15.“ Thrust and power characteristics for constant forward velocity; 
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Figure 15.- Concluded, 
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